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Recently it has been shown that time optimal quantum computation is attained using the Cartan 
decomposition of a unitary matrix. We extend this approach by noting that the unitary group is 
compact. This allows us to reduce the execution time of a quantum algorithm U a \ g further by adding 
an extra gate W to it. This gate W sends ?7 a i g to another algorithm WUa\ s which is executable in 
a shorter time than t/ a i g . We call this technique warp-drive. Here we show both theoretically and 
experimentally that the warp-drive reduces the execution time of Grover's algorithm implemented 
with a two-qubit NMR quantum computer. Warp-drive is potentially a powerful tool in accelerating 
algorithms and reducing errors in any realization of a quantum computer. 

PACS numbers: 03.67.Lx, 82.56.Jn 

Keywords: warp-drive quantum computation, time optimal control, Cartan decomposition, Grover's algo- 
rithm, NMR quantum computer 



Quantum computing is an emerging discipline based on 
encoding information into a quantum-mechanical system 
P, ■ There a quantum algorithm is expressed in a form 
of a unitary matrix 



E/aig = Texp 



H(-y(t))dt 



(1) 



where T stands for time-ordered product, 7(f) denotes 
collectively the control parameters of the Hamiltonian 7i 
at time t. This unitary matrix is often implemented in 
terms of so-called elementary gates, such as U(2) gates 
and CNOT gates [J It is possible, instead, to 
directly implement a given unitary matrix without de- 
composing it into these elementary gates. It is expected 
that this will reduce the execution time required in gen- 
eral. One method for direct implementation is to employ 
numerical optimization of the control parameters of the 
Hamiltonian 0, Q> @ • The other method is to use the 
Cartan decomposition of the group SU(2 n ), to which an 
7i-qubit matrix representation of a quantum algorithm 
belongs [jj ^| . The latter approach, which is adopted 
here, is successfully demonstrated recently using an NMR 
quantum computer, whose pseudopure state is generated 
by cyclic permutations of state population [TlT | . It should 
be noted that exact optimal implementation of a quan- 
tum algorithm has been achieved in holonomic quantum 
computation in an idealized case |12| . 

A warp-drive is a fictitious gadget with which two re- 
mote points in space are connected |l3l . It is the purpose 
of the present Letter to demonstrate, both theoretically 
and experimentally, that a similar technique may be em- 
ployed to shorten the execution time of a quantum algo- 
rithm. A time optimal implementation of the quantum 



algorithm is equivalent to navigating along the time op- 
timal path from the identity operator to the point U a \ g 
in SU(2") by tuning the parameters -f(t) in the Hamil- 
tonian. We show below that an additional permutation 
matrix W of the basis vectors, when added after U a \ gl 
sends it to a point WU a \ g near the identity matrix / 
so that it takes a shorter time to follow the time opti- 
mal path connecting I with WU a \ g than / with J7 a i g , see 
Fig. 1. Therefore we have "warp-driven" U a \ g to WU a \ g 
by adding W . Although it might seem counterintuitive 
that the execution time is reduced by adding an extra 
gate, this is the case since the unitary group SU(2") is 
compact. Let us consider navigating on a sphere S 2 , 
which is also a compact space. Suppose we leave from 




FIG. 1: Conceptual diagram showing the effect of a warp- 
drive gate W in a compact space SU(4). A matrix [7 a i g is 
sent to WUaig which is reachable from the unit matrix 7 in a 
shorter execution time. The curves connecting these matrices 
represent time optimal paths. 
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a point on the equator toward west. Then we will 
eventually arrive at the antipodal point. The distance 
between O and us will be shorter and shorter as we fur- 
ther circumnavigate the sphere. This is what happens 
when an extra gate is added. 

In an NMR quantum computer, a one-qubit operation 
may be carried out in a short time on the order 10 /its 
while a two-qubit entangling operation takes time typ- 
ically ~ 10 ms. Thus one-qubit operation time may be 
ignored in estimating the overall execution time. Let us 
consider a molecule with two heteronucleus spins for def- 
initeness, whose Hamiltonian, in the rotating frame with 
respective Larmor frequency, is 

~H{-f) = — wn [cos fa (<x x ® 72/2) + sin fa (a y <g> J2/2)] 
— W12 [cos fa(h ® <y x /2) + siYL(j) 2 {I 2 ® 0^/2)] 
+27rJcr 2 ® er 2 /4, (2) 

where wij and ^ are the control parameters denoting 
the amplitude and the angle in the xy-plane of the exter- 
nal rf fields, respectively. Above J denotes the spin-spin 
coupling constant. Typically we have uiu 3> J, which 
justifies the above assumption of negligible one-qubit op- 
eration time compared to two-qubit operation time. Ne- 
glecting one-qubit operation time in evaluating the exe- 
cution time amounts to identifying the matrices U± and 
U 2 which differ by an element of K = SU(2) <g> SU(2). 
Thus the relevant space for evaluating the time opti- 
mal path is the coset space SU(4)/SU(2) <g> SU(2). To 
find the time optimal path connecting the unit matrix I 
and the matrix C/ a i g , therefore, amounts to finding the 
time optimal path connecting cosets [I] and [£7 a i g ], where 
[U] = {kU\k e K}. The Lie algebra su(4) of SU(4) is 
decomposed as su(4) = t © p 0] , where 

6 = Span({iJ(g)(7j/2,ia i (g)I/2}), (j = x,y,z) (3) 
P = 6 X = Span({icr J - ® cr fc /4}), (J, k = x, y, z). (4) 

Note that they satisfy the commutation relations 

[t,f]ct, [p,t]cp, [p,p]ce. 



(5) 



The decomposition of a Lie algebra g into t and p, sat- 
isfying the above commutation relation, is called the 
Cartan decomposition. The Cartan subalgebra f) = 
Span({i<7j £g> <7j/4}) C p plays an important role in the 
following construction. A general theorem of Lie algebras 
proves that any element U a i g € SU(4) has a KP decom- 
position C/ a ig = kp, k € K = exp6 and p G P = expp. 
Moreover, any matrix p G P is rewritten in the conjugate 
form p — k\hk\, where k\ G K and h is an element of 
the Cartan subgroup H of SU(4), 



H = cxp f) = < cxp I i ^ ~A~ a i ® a i 



(6) 

Therefore we have a corresponding Cartan decomposi- 
tion for a group element as [/ a ig = k 2 hki, where fcj G K 



and h E H. The implementation of a quantum algo- 
rithm based on a Cartan decom pos ition requires shorter 
execution time in general 0, ITflll l| . 

To be more concrete, let us consider implementing two- 
qubit Grover's database search algorithm JTy with an 
NMR quantum computer. The data is encoded in one 
of the basis vectors |00), |01), 1 10) , |11) and the gate Uij 
picks out a p articular binary basis vector \ij) as a "target 
file" hj,hj. Here we restrict ourselves within Uio which 
picks out the file 1 10) with a single step. The unitary 
matrix representing this algorithm takes the form 



Ui = 






1 



















-1 

















-1 






(7) 



The Cartan decomposition of an arbitrary U G SU(4) is 
carried out explicitly as follows. We first introduce the 
Bell basis E3 



|* ) = (1A/2)(|00) + |11)), 

|*!> = (</>/S)(|oi) + |io», 

|* 2 ) = (1/V2)(|01)-|10», 

|* 3 ) = (VV2)(|oo)-|H)). 



(8) 



The transformation rule of a matrix U with respect to the 
standard binary basis |00), 1 01) , |10), 1 1 1) into that with 
the Bell basis is U -> U B = Q^UQ, where 



1 
71 



( 1 i 

i 1 

i -1 

V 1 -i 



(9) 



The matrix Q defines an isomorphism between K = 
SU(2) <g> SU(2_^ and SO(4) and is used to classify two- 
qubit gates [TtI Il8| . Namely, it is easy to verify that 
Q^kQ G SO(4) for k G K. Moreover, Q diagonal- 
izes the elements of the Cartan subgroup, viz Q^hQ = 
diag(e te °,e w \e iS2 ,e W:3 ) for h G H. Therefore we find for 
U = k 2 hki that 

U B = Q^UQ = Q^k 2 Q ■ Q^hQ ■ Q^k 2 Q = 2 h D O l , 

where Oi = Q^kiQ G SO(4) and ho = Q^hQ is a diag- 
onal matrix. From U^Ub = Ofh D Oi, we notice that 
UbIIb is diagonalized by 0\ and its eigenvalues form 
the diagonal elements of h D . Finally 2 is found as 

o 2 = u B {h D o l y 1 . 

We apply the above strategy to find the Cartan decom- 
position Uio = k 2 hk\. An example of the time optimal 
control is 



h ® h, 

i(ir/4) (<r x (g>rr x -<7 y ®o-j, ) 



ki 
h 



(10) 



To implement this decomposition with an NMR quantum 
computer, such terms as e 2 ^/ 4 ^ '*® ' 1 ) must be rewritten 
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TABLE I: time optimal pulse sequences for Grover's algorithm Uio and the warp-driven algorithm, W4U10. The hydrogen 
nucleus is the qubit 1 while the carbon nucleus is the qubit 2. Here X (Xm) and Y (Ym) denote 7r/2-pulse around x (— x) and 
y (—y) axis, respectively. The symbol Pi(6) denotes a 7r-pulse around a vector (cos 0, sin 0, 0) in the Bloch sphere. The symbol 
(1/2 J) indicates the length of the idle time, during which no external pulses are applied. The number of pulses is reduced from 
10 to 4 and the execution time is halved by adding the extra gate Wa. 



Gate 




Pulse sequence 


Execution Time 


t/10 


1: 


X (1/2J) Xm Y (1/2J) X 


Ym 1/J 




2: 


X (1/2J) Xm Ym (1/2J) Y 


PiO/4) 


W 4 U W 


1: 


Xm Pi(-7r/4) (1/2J) X 


1/2J 




2: 


(1/2J) Pi(77) 





in favor of the subset of generators of SU(4) in the Hamil- 
tonian (1). We verify, for example, that 



r e i(iv/i)a x e -i(ir/4)a y -i e i(ir/4)(<Tz<g>cT z ) 



(11) 



The extra gate W which possibly shortens the execu- 
tion time must be simple enough so that we can deduce 
the output of Ugig efficiently from that of WU a i g using 
classical computation only. We call such gates warp-drive 
gates. Since a matrix U a i g is an element of a compact 
group U(2"), there always exist such warp-drive gates 
which will reduce the execution time. In the present 
work, we choose the permutation matrices of the binary 
basis vectors as candidates for such gates. Since there 
are four basis vectors, there are 4! = 24 permutation ma- 
trices. Note that we are not required to examine all of 24 
permutations since 18 of them are obtained by applying 
one-qubit rotations on the following six permutations: 



W 
W 5 



^SWAP, W4 — U cp = C^CNOT^CNOT' (12) 



TJ 21 TJ 12 

u CNOT u CNOT J 



CNOT 



IS 



where J7swap represents the swap gate and U, 
the CNOT gate whose control bit is i while the target 
bit is j. The matrices U cp and U 2 p have been utilized 
to cyclically permute the state populations to generate 
pseudopure states [l9|. 

We have optimized the execution time of the matrices 
WiUio (0 < i < 5) by utilizing the Cartan decomposition 
outlined above and found that the execution time is 1/ J 
for i = 0, 1 and 2 while 1/2 J for i — 3,4 and 5. A time 
optimal control for the warp-driven gate W4U10 = kihk\ 



h = e i ^/ A ) a v® a v 



(13) 



ko = If 



TCx/4 



Therefore the execution time T satisfies irJT/2 = tt/4, 
yielding T = 1/2 J. Since the decomposition l(T^|) con- 
tains generators which do not exist in the Hamiltonian, 
a trick similar to (|11|) must be employed. The one-qubit 
gates are realized using rf-pulses in an NMR quantum 



computer. Table I shows the actual NMR pulse se- 
quences derived from the Cartan decompositions of U\q 
and W4U1Q. We call the hydrogen nucleus and the car- 
bon nucleus as qubit 1 and qubit 2, respectively, in Table 
I and the rest of this Letter. In spite of the extra gate 
W4, the latter requires less pulses and half of the execu- 
tion time required for the former. The permutation W4 
maps the binary basis vectors as W4 : 1 10) 1— * |11), and 
accordingly the output state of W4II10 upon acting the 
initial state |00) is |11). 

In our experiments, a 0.6 ml, 200 mM sample of 
carbon- 13 labeled chloroform (Cambridge Isotope) in d-6 
acetone has been employed as a two-qubit molecule and 
data were taken at room temperature with a JEOL ECA- 
500 NMR spectrometer, whose hydrogen Larmor fre- 
quency is approximately 500 MHz. The measured spin- 
spin coupling constant is J = 215.5 Hz and the transverse 
relaxation time is T2 ~ 7.5 s for the hydrogen nucleus and 
T2 ~ 0.30 s for the carbon nucleus. The longitudinal re- 
laxation time is measured to be T\ ~ 20 s for both nuclei. 
The initial state |00) is prepared as a pseudopure state 
generated by the field gradient method poj . 

Figure 2 shows our experimental result. The spectra 
for the 13 C nucleus (qubit 2) are shown in the panels. The 
spin states of both nuclei are inferred from these spectra. 
The position of the peak depends on the state of the qubit 



1 00) 



|10) 



tyio> - 



' 79.20 ppm 



77.49 ppm 
11) 



FIG. 2: NMR spectra obtained by applying a reading pulse to 
the 13 C nucleus. The rear row shows the pseudopure initial 
state |00) while the middle row indicates the state 1 10) ob- 
tained by executing U\o on 1 00) . Finally the front row shows 
the spectrum obtained from warp-drive quantum computing 
WaUio, where the peak corresponds to |11). 
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1: peaks at 77.49 ppm (part per million) and 79.20 ppm 
stand for the states |1) and |0), respectively, of the qubit 
1. The positive (negative) amplitude indicates qubit 2 
being in the state |0) (|1)). Thus the peak in the rear row 
shows the initial pseudopure state is |00) while the peak 
in the middle row tells us that the execution of Grover's 
algorithm U±o generates the state |10). Finally the front 
row shows the spectrum obtained after the execution of 
the warp-driven gate W4U10. The peak shows that the 
resulting state is |11), from which the output of J7io|00) is 
easily deduced as |10). The spectrum obtained with the 
warp-driven gate is sharper than that with the original 
Uio gate. The undesirable peak at |00) is also improved 
by the warp-drive. Smaller number of gates and reduced 
execution time account for these improvements. 

We have also examined the effect of the warp-drives 
for Grover's algorithms other than Uio and found that 
the execution times remain unchanged for Wo, 1,2 while 
they are halved for ^3,4,5. It should not be expected, 
however, that the latter gates reduce the execution time 
of an arbitrary two-qubit gate: clearly, such a "universal 
warp-drive" which reduces execution time of an arbitrary 
quantum algorithm does not exist. In view of the com- 
pactness of the unitary group, however, there exists a fi- 
nite set of simple unitary transformations, such that the 
execution time of any quantum algorithm is reduced by 



properly chosen elements of the set. We expect that the 
set of permutation matrices, or certain subset thereof, is 
such a "universal warp-drive set" . 

In summary, we have demonstrated both theoretically 
and experimentally that a quantum algorithm may be ac- 
celerated by adding an extra gate to it. We took advan- 
tage of the compactness of the group SU(4) and found the 
time optimal control parameters utilizing the Cartan de- 
composition for it. When applied to two-qubit Grover's 
algorithm in the NMR experiment the execution time is 
found to be halved along with a reduction in the number 
of the gate pulses from 10 to 4. As a result, application of 
warp-drive quantum computing sharpens the NMR spec- 
trum and reduces the spurious peak, implying reduction 
in decoherence and gate operation errors. How reduc- 
tions in the execution time and the number of gate pulses 
scales with the number of qubits remains a challenging 
problem. 
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